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We calculate the dynamical structural factor of the S = 1 bond-alternating Heisenberg chain. In 
the Haldane phase, the lowest excited states form the lower edge of the multimagnon continuum in 
< 9 < gc and the one-magnon mode 'm. qc < q < n. As the system approaches the gapless point, 
qc shifts towards q = n and the largest integrated intensity of the one-magnon mode is decreased. 
In the singlet-dimer phase, the one-magnon mode appears in < g < Qc- As the bond-alternation 
becomes strong, qc shifts towards q = n. In the antiferromagnetic-ferromagnetic bond-alternation 
region with a strong ferromagnetic coupling, the lowest excited states form the lower edge of the 
multimagnon continuum in < 5 < 0.27r and O.Stt < </ < tt, and the one-magnon mode appears in 
0.27r < q < O.Stt. The largest integrated intensity of the one-magnon mode is 93%, which is slightly 
smaller than that in the 5 = 1 Haldane-gap system. We further discuss the dynamical structural 
factor in connection with the inelastic neutron-scattering experiments. 

PACS numbers: 75.40.Gb, 75.10.-b, 75.40.Mg 



I. INTRODUCTION 

The isotropic Heisenberg chain with integer spins has 
a finite excitation gap above a disordered ground state, 
while that with half-odd-integer spins has a gapless ex- 
citation from a critical ground state 0, This Hal- 
dane conjecture based on the 0(3) nonlinear a model 
with a topological term has been confirmed by exten- 
sive numerical and experimental studies. The elemen- 
tary excitations of the S = 1 isotropic Heisenberg chain 
were investigated by inelastic neutron-scattering exper- 
iments. Using S — 1 quasi-one-dimensional Heisen- 
berg antiferromagnets Ni(C2H8N2)2N02C104 and 
CsNiCla 5j Ji], the dynamical structure factors S{q,LLj) 
for the one-magnon mode and the multimagnon contin- 
uum were observed. The integrated S{q, uj) with respect 
to a; at g = TT from the multimagnon mode was compared 
with the results obtained by several theoretical methods 
nil, H EIH m ■ The origin of the deviation was dis- 
cussed in connection with the picture for the elementary 
excitation Q. 

After some years later from the Haldane conjecture, 
Affleck [H and Affleck-Haldane ^ further applied 
the same analysis to the bond-alternating spin-S* Heisen- 
berg chain described by the following Hamiltonian: 

H = J ^ 1^82,-1 • S2. + aS2. • S2,+i^ , (1) 

and argued that there exist 25' gapless points in < a < 
00. Their argument was successfully confirmed by several 
theoretical methods [IlllllllllillSIIlliiillll. 

For the S = 1 system, the critical a separating the 
Haldane phase {ac < a < 1) and the singlet-dimer 
phase (0 < a < a^) was estimated using a series ex- 
pansion method |l6ll. a density-matrix renormalization- 
group method |l7j. quantum Monte Carlo methods 0, 
Uyj I23L , the Binder parameter combined with finite- 
size scaling and a level-crossing method l25|. The 



results are in agreement as ac ~ 0.6. When — cx3 < a < 
in the S = 1 system, on the other hand, the trivial dimer 
phase around a ~ is shown to be smoothly connected 
to the S = 2 Haldane phase which appears in the limit 
a -> -00 [23. 

Experimentally, using some nickel chain compounds 
which are regarded as the S = 1 bond-alternating Heisen- 
berg system, magnetic susceptibility, ma gnet ization, and 
ESR measurements were performed jH 11 El El El . 
From the comparison between the experimental and nu- 
merical results, the strength of the bond-alternation a 
was evaluated and it was concluded that some com- 
pounds are in the singlet-dimer phase psl l26 l. '28'| and 
the other compound is close to the gapless point 23, 27j. 

The elementary excitations of the S = 1 bond- 
alternating Heisenberg chain were investigated by a 
quantum Monte Carlo method and an exact diag- 
onalization method 20]. From the dispersion relation 
of the low-lying excitation, it was shown that the gap 
opens at the center of the Brillouin zone for a > 0, while 
the gap opens at the boundary of the Brillouin zone for 
a < 0. In the Haldane phase around a = 1, the low- lying 
excitations are expected to be scattering states of the do- 
main walls in the hidden antiferromagnetic ordering |19| . 
Around a = 0, in contrast, the low-lying excitations are 
well described by the S — 1 magnon j23|. In spite of 
these studies, dynamical properties have not yet been 
fully investigated. A detailed investigation is desirable. 

In this paper, we investigate dynamical properties of 
the S = I bond-alternating Heisenberg chain. Using 
a continued fraction method based on the Lanczos al- 
gorithm El, we systematically calculate the dynamical 
structure factor (DSF) in the Haldane phase and dimer 
phase including the regions a > and a < 0. The DSF 
provides us with information about the intensity of the 
magnetic excitation as a function of energy-momentum 
transfer. Since the DSF can be observed by inelastic 
neutron-scattering experiments, a comparison between 
theoretical and experimental results is possible. The 
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setup of the paper is as follows. In Sec. II, we briefly 
summarize the method for the numerical calculation. In 
Sec. Ill, we show the results for the distribution of the in- 
tensity as a function of energy- momentum transfer. Us- 
ing finite-size effects @, we discuss whether the lowest 
excited states in given wave numbers form an isolated 
mode or a lower edge of the excitation continuum. Sec. 
IV is devoted to the summary of the paper. 



II. MODEL AND METHOD 

We consider the = 1 bond-alternating Heisenberg 
chain described by the Hamiltonian (1). We set that 
the total number of spins is A'^ = 20 at the maximum 
and J — 1. We apply the periodic boundary condition. 
The DSF can be expressed in the form of the continued 
fraction as 

TT ^ z — H ^ 

= S^{q)C'^{q,cu) {fi^x,y,z), (2) 

where Iv^o) is the eigenstate with the lowest eigenvalue 
Eo, = {llVN)Y.je""S'^ and z = tu + ir] + Eq 
with h = 1. The lattice constant between neighboring 
two sites is set to unity. Therefore, q = 0.5tt corre- 
sponds to the boundary of the Brillouin zone. In the 
expression (2), S'^{q) is the static structure factor and 
C^{q,u}) is represented in the form of the continued 
fraction, which can be calculated numerically by Lanc- 
zos algorithm. The total contribution of C'^{q,uj) for 
fixed q is normalized to unity, because the following sum 
rule has to be satisfied: S'^{q) = duiS'^{q,Lo). Since 
rotational symmetry around a;, y and z axes remains, 
S^{q,U!) = S^iq^Lo) — S^{q,uj). We thus calculate only 
S^{q,cj) = S{q,Lj){S^{q)C^{q,iu) = S{q)C{q,u;)). In- 
stead of taking 77 +0, we set rj = 1.0 x 10 ^. For 
finite-size systems, therefore, C^{q,uj) consists of a finite 
number of Lorentzians. Following Ref. 9, we call the 
position of a Lorentzian and the integrated value of each 
Lorentzian with respect to u as the pole and residue, 
respectively. 

To discuss whether the lowest excited states in given q 
form an isolated mode or a lower edge of the excitation 
continuum in the thermodynamic limit, we investigate 
the finite-size effects of the poles and their residues of 
the continued fraction C{q,Lo) 0,1131 . As reported in 
Refs. 9 and 30-32, a pole which belongs to an excitation 
continuum tends to have an appreciable size dependence 
on at least either its position or its residue, while a pole 
which belongs to an isolated mode hardly depends on the 
system size. 

Using these methods, we proceed to the numerical cal- 
culation for the DSF and discuss characteristics of the 
low-lying excitation. 



III. RESULTS 

We first investigate characteristics of the lowest excited 
state in a given q. In Fig. 1, the finite-size effects for the 
residues of the lowest excited states are shown in the ex- 
tended zone scheme. Note that the finite-size effects only 
for the residues are presented, since the positions of the 
poles scarcely depend on the system size. We discuss the 
results, making reference to the typical size dependence 
of the excitation continua for the 5=1/2 's^ and 1 
isotropic Heisenberg chains. 

In Fig. 2, the DSF's of the Haldane phase and the 
dimer phase are shown in the extended zone scheme. The 
intensity of S^{q,uj) is proportional to the area of the 
circle. The full circles represent the isolated mode and 
the gray circles represent the excitation continuum. The 
largest intensity in a given q appears at the lowest excited 
state. In a > it increases as q approaches tt, while in 
a < it increases as q approaches O.Stt. 



A. Haldane phase 

As N increases, the residues at a = 0.8 and 0.7 de- 
crease in < g < O.Stt and in < g < O.Stt, respec- 
tively. The size dependence of the residues in both re- 
gions of wave numbers is quite similar to that for the ex- 
citation continua of the 5 = 1/2 jSfl] and 1 ^ isotropic 
Heisenberg chains. The behavior of the residue around 
q = 0.9tt at a = 0.7 is different from that in the other 
wave numbers. Such size dependence is probably caused 
by the fact that the next lowest excited state around 
q = 0.9tt lies close to the lowest excited state with almost 
the same scattering intensity. Therefore, the size depen- 
dence around q = 0.9tt suggests that the lowest excited 
states may become the lower edge of the multimagnon 
continuum. Accordingly, at a = 0.8, the lowest excited 
states in < g < O.Stt become the lower edge of the mul- 
timagnon continuum, while those in 0.5tt < q < it form 
the one-magnon isolated mode. At a = 0.7, the lowest 
excited states in < 5 < 0.97T become the lower edge of 
the multimagnon continuum, while those around q = tt 
form the one-magnon isolated mode. 

In Fig. 2, we show the DSF at a = 0.8 and 0.7 in 
the extended zone scheme. The largest intensity caused 
by the one-magnon mode appears at g = tt. We evalu- 
ate its integrated intensity by extrapolating at A'^ — > 00. 
The ratio that the one-magnon mode carries at q = tt is 
87% in a = 0.8 and 74% in a = 0.7, respectively. These 
values are smaller than that in the 5=1 Haldane-gap 
system, 97% ~ 98% Hill [Hill. Therefore, as the sys- 
tem approaches the gapless point in the Haldane phase, 
the one-magnon mode is reduced towards g = tt with de- 
crease in the excitation gap and the scattering intensity 
at g = TT. 

At a = 0.6. the system is in the vicinity of the gapless 
point [SllllllliliSlllllllll. The residues in 
< g < TT decrease with increasing iV, indicating that 
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the lowest excited states in < g < tt form the lower 
edge of the excitation continuum. 



Relation to inelastic neutron-scattering 
experiments 



B. Dimer phase 

Judging from the finite-size effects of the residues, the 
one-magnon isolated mode appears in < q < O.Stt at 
a — 0.5, while the lowest excited states form the one- 
magnon mode inO<(7<7ratQ; = 0.3 and 0.1. As 
the system leaves from the gapless point in the singlet- 
dimer phase, the wave-number region of the one-magnon 
mode seems to extend towards q = n. At a = —0.3, 
the lowest excited states form the one-magnon mode in 
< g < O.Stt and the lower edge of the multimagnon 
continuum in O.Stt < g < tt. 

At a = —10.0, the system can be regarded as the 5 = 2 
Haldane-gap system. From the finite-size effects of the 
residues, we find out that the one-magnon mode appears 
in 0.27r < g < O.Stt, and the lowest excited states form 
the multimagnon continuum in < g < 0.2tt and O.Stt < 
g < TT. The integrated intensity of the one-magnon mode 
at g = O.Stt is evaluated to be 93% by extrapolating at 
N oo. The integrated intensity that the one-magnon 
mode carries at the boundary of the Brillouin zone in the 
S = 2 Haldane-gap systems is smaller than that in the 
S = 1 Haldane-gap system. 
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FIG. 1: The finite-size effects of the residues of the lowest 
excited states in given q. 



We now investigate the observable DSF by inelastic 
neutron-scattering experiments. The excitation energies 
of the lowest excited states are summarized in Table L 
Comparing the excitation energies at the symmetric wave 
numbers about g = 0.5tt, we discuss the DSF in the re- 
duced zone scheme. 

At a = O.S and 0.7 the excitation energies in O.Stt < 
g < TT lie below the corresponding ones in < g < O.Stt. 
Therefore, at a = O.S, the lowest excited states in 
< g < O.Stt in the reduced zone scheme form the one- 
magnon isolated mode. At a = 0.7, the lowest excited 
states around g = form the one-magnon isolated mode, 
while those in O.Itt < g < O.Stt form the lower edge of the 
multimagnon continuum. The largest intensity caused 
by the one-magnon mode appears at g = with the inte- 
grated intensity ~ 87% in a = O.S and ~ 74% in a = 0.7. 

In the dimer phase, the excitation energies take the 
same values at the symmetric wave numbers about g = 
O.Stt within our numerical accuracy. Thus, at a = 0.3 
and 0.1, the lowest excited states in the reduced zone 
scheme form the one-magnon mode in < g < O.Stt. 
At a — O.S and —0.3, the one-magnon mode coincides 
with the lower edge of the multimagnon continuum. In 
this case, the one-magnon mode may become unsta- 
ble and only the multimagnon continuum emerges. In 
fact, it was observed by inelastic neutron-scattering ex- 
periments that, when the one-magnon mode crosses the 
multimagnon continuum, the large intensity of the one- 
magnon mode becomes invisible in the multimagnon con- 
tinuum [1,0,IE@1- At a = O.S, therefore, the lowest 
excited states form the lower edge of the multimagnon 
continuum in < g < O.Stt. At a = —0.3, the lowest 
excited states form the lower edge of the multimagnon 
continuum in < g < 0.2tt and the one-magnon mode in 
0.27T < g < O.Stt. 

At a = 0.3 and 0.1, the one-magnon modes corre- 
sponding to the symmetric wave numbers coincide with 
each other. At a = —0.3, there exists such region of the 
wave numbers. In these cases, the intensity is described 
as the sum of two one-magnon modes. Therefore, the 
largest intensity does not always appear at g = or O.Stt 
in the reduced zone scheme. At a = 0.3 and 0.1, the in- 
tensity of the one-magnon mode increases as g approaches 
0. At a = —0.3, in contrast, the largest intensity of the 
one-magnon mode emerges at g = 0.4tt. 

At a — —10.0, the multimagnon continuum appears 
in < g < 0.2tt and the one-magnon mode appears in 
0.2Tr < g < O.Stt with the largest scattering intensity at 
g = O.Stt. The wave-number region for the one-magnon 
mode is smaller than that in the S = I Haldane gap sys- 
tem Q . The largest intensity caused by the one-magnon 
mode appears at g = O.Stt with the integrated intensity 
93%. 

The DSF's in the reduced zone scheme are presented 
in Fig. 3. Experimentally, thermodynamic properties of 
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the S = I bond-alternating Heisenberg chain have been 
investigated so far. To obtain detailed information on 
the elementary excitation, the study of dynamical prop- 
erties is needed. Such studies should shed more light on 
understanding the nature of the S = 1 bond-alternating 
Heisenberg chain. 
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FIG. 2: S{q.,u!) for A'^ = 20 in the extended zone scheme. 
The intensity is proportional to the area of the circle. The 
full circles represent the isolated branch and the gray circles 
represent the excitation continuum. 



< 9 < (Zc- As the bond- alternation becomes strong, qc 
shifts towards q = n. At a = —10.0, the lowest excited 
states form the lower edge of the multimagnon continuum 
in < g < 0.2-77 and O.Stt < q < tt, and the one-magnon 
mode appears in 0.2^ < q < O.Stt. 

We have further discussed the DSF in the reduced 
zone scheme in connection with the inelastic neutron- 
scattering experiments. 
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FIG. 3: S(q.,uj) for A?^ = 20 in the reduced zone scheme. 
The intensity is proportional to the area of the circle. The 
full circles represent the isolated branch and the gray circles 
represent the excitation continuum. 



IV. SUMMARY 
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TABLE I: The excitation energies of the lowest excited states 
in the extended zone scheme for N = 20. The boundary of 
the Brillouin zone is g = O.Stt. 

Q q = n/10 27r/10 Stt/IO 47r/10 Stt/IO 6n/10 Ttt/IO Stt/IO Qtt/IO IOtt/IO 
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